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ABSTRACT 

We present the power spectrum of galaxy clusters measured from the new ROSAT- 
ESO Flux-Limited X-Ray (REFLEX II) galaxy cluster catalogue. This new sample 
extends the flux limit of the original REFLEX to 1.8 x 10 -12 ergs -1 cm -2 , yielding a 
total of 911 clusters with ^ 94 per cent completeness in redshift follow-up. The analysis 
of the data is improved by creating a set of 100 REFLEX Il-like mock galaxy cluster 
catalogues built from a suite of large volume ACDM iV-body simulations (L-BASICC 
II). The measured power spectrum is in agreement with the predictions from a ACDM 
cosmological model. The measurements show the expected increase in the amplitude 
of the power spectrum with increasing X-ray luminosity. On large scales, we show that 
the shape of the measured power spectrum is compatible with a scale independent bias 
and provide a model for the amplitude that allows us to connect our measurements 
with a cosmological model. By implementing a luminosity-dependent power spectrum 
estimator, we observe that the power spectrum measured from the REFLEX II sample 
is weakly affected by flux-selection effects. The shape of the measured power spectrum 
is compatible with a featureless power spectrum on scales k > 0.01 h Mpc -1 and hence 
no statistically significant signal of baryonic acoustic oscillations can be detected. We 
show that the measured REFLEX II power spectrum displays signatures of non-linear 
evolution. 
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1 INTRODUCTION 

Within the last decade, three foundational observational 
probes have been well recognized as opening up observa- 
tional windows to reveal some of the most valuable se- 
crets of the Universe. These are: the temperature fluc- 
tuations of the cosmic microwave background radiation 
(CMB), recently measured with high precision by the 
Wilkinson Microwave Anisotropy Probe ( WMAP) satellite 
(|Spergel et all 120071 : iKomatsu et all 120101 ); the Hubble di- 
agram inferr ed from the Type la supernovae (SNe) o bser- 
vations (e.g. IPerlmutter et all 1 19991 : iRiess et al.l 120041 ) and 
the measurement of the large scale structure (LSS) of the 
Universe as traced by the spatial distribution of galaxies 



dPercival et all 12002]; iTegmark et al. I 120061; ISanchez et all 
l2006tlPercival et alj|2007l l20ld : iReid et aUl201Ch . 

With the completion of large redshift surveys, such as 
the Two-degree Field Galaxy Redshift Survey (2dFGRSfl 
and the Sloan Digital Sky Survey (SDSSfl it has been pos- 
sible to push the level of accuracy of LSS studies. The re- 
cen t detection of the baryon acoustic oscillations (BAO, 
e.g. Eisens tein et al.ll2005l:ICole et a l. 2005; Gaztanag a et al.l 



l200St ISanchez et alj|2009l : iPercival et alj l201(t > in redshift 
surveys have been key to this progress. Over the last few 
years these observations have established the concordance 
cosmological model, based on a flat space-time in a current 
phase of accelerated expansion due to the presence of a dom- 
inating dark energy component, whose equation of state is 
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compatible with Einstein's cosmological constant A. This is 
the so-called ACDM cosmological model. 

Galaxy cluster samples also provide valuable infor- 
mation to obtain constraints on cosmological parameters. 
Galaxy clusters are the largest bounded structures in the 
Universe. They are associated with the highest peaks in the 
matter density field and are recognized as biased tracers 
of th e underlying matter distribution (e.g. iBardeen et al.l 
1986). Their deep potential wells make them the largest 
astrophysical laboratories in the Universe, where the com- 
bination of gravitation and baryonic physics has been in- 
tensively studied through the analy sis of cluster proper- 
ties s uch as scaling relations (e.g. Reiprich fe Bohringer] 



I199E 



Pratt et al 



Makino et all 



20091 ; iMantz et al.ll2010t ), d ensity profiles 



1998 ) , pres sure profiles (e . g.lArnaud et all 
l2009h . baryon fractions (e.g. iGiodini et afll2009h . etc. The 



abundances of ga laxy clusters determin ed by their luminos- 
ity function (e.g. iBohringer et ;il.ll2llll2 :■ can also be used to 
constrain parameters like the matter content in the universe 
ttm, and the amplitude of the density fluctuations charac- 
terized by 08, the rms linear pertur bation theory varianc e 
in spheres of radius 8 Mpc h~ l (e.g. ISchuecker et al.| [2003). 
The spatial distribution of galaxy clusters, characterized by 
its power spectrum (or correlation function), provides useful 
information about the cosmological model of the Universe. 
The shape of this measurement is particularly sensitive to 
the parameter combination fl m h, which complements the 
constraints on these parameters obtained from the analysis 
of fluctuations in the CMB. Furthermore, the amplitude of 
the galaxy cluster power spectrum also contains important 
information that can be related to theoretical models in a 
more direct way than for measurem ents based on galaxy 
samples (e.g. [Moscardini et al. 200(J). 

In the past few years t he ROSAT-ESO Flux L imited X- 
ray (REFLEX) catalogue (|Bohringer et al.ll200ll ) has been 
used to measure fundamental cosmological quantities. The 
REFLEX catalogue i s based on ROS AT All Sky Survey 
(RASS) observations (|Truemper Ill993 |) , complemented wit h 
follow-up observations as described bv lGuzzo etafl l|2009l) . 
yielding spectroscopic redshifts for 447 clusters with flux 
limit of 3 x 10~ 12 ergs _1 cm~ 2 (in the ROSAT energy band 
0.1 - 2.4 keV). The REFLEX catalogue was, to date, the 
largest statistically complete X-ray detected cluster sam- 
ple. The clustering properties of thi s survey were analyze d 
by means of the power spectr um (ISchuecker et al] l200ll ) , 
the cluster correlation function (ICollins et al. Il2000h. cluster- 
galaxy cross-correlation functions { Sanchez et al. 2005]) and 
Minkowski functionals (|Kerscher et al.l|200ll ). Sub-samples 
of the REFLEX catalogue complemented by detailed follow- 
up observations have been used to constrain cluster scal- 
ing relations (e.g. iReiprich fe Bohringerl|l999l ; IStanek et al.l 
120061 ; IPratt et al.ll2009l ; IMantz et al.ll2010l ). 

In this paper we present the analysis of the power spec- 
trum of the new REFLEX II catalogue. The REFLEX II is 
an extension of the REFLEX catalogue to a lower limiting 
flux (1.8 x 10 -12 ergs -1 cm" 2 ) allowing the inclusion of 464 
new clusters over the original sample and yielding a total of 
911 clusters with spectroscopic redshifts for ~ 95 per cent 
of the sample. 

In addition to the enlarged sample size of REFLEX II, 
several improvements were made to the data reduction: (i) 
we use the RASS survey product RASS III which gives a 



few percent more sky exposure in formerly underexposed 
areas due to improved attitude solutions, consequently re- 
covering a few more clusters at higher flux, (ii) for the count 
rate to flux conversion an estimated temperature has to be 
applied, which is now obtained with up-to-date scaling rela - 
tions based on the REXCESS Surve y llBohringer et al|2007l) 
with the L-T relation described in IPratt et al.l 1)20091 ); (hi) 



the total flux and X-ray luminosity is now estimated in- 
si de the radii of rsq o an d r2oo (based on rel ations described 
in lPratt et aill2009l and lArnaud et al.ll2005l ). These calcula- 
tions now involve less extrapolation than the estimates for 
the previously used fiducial radius. 

Besides the advantages provided by a larger cluster 
sample, the power spectrum analys is presented here repre- 
sents an improvement over that of ISchuecker et all (|200ll ) 
in a number of ways. In particular, our analysis is comple- 
me nted with a set of A f-body simulations, t he L-BASICC 
II (|Angulo et all 120081 ; ISanchez et "all l2008al ). from which 
we constructed a suit of 100 REFLEX II mock catalogues. 
These catalogues were calibrated to reproduce the measured 
REFLEX II X-ray luminosity function. Selection criteria of 
the REFLEX II sample were applied in their construction, 
yielding a large suit of mocks that can be used to analyze 
the statistical methods applied to the data. The details of 
the construction of these mock catalogues will be described 
in a forthcoming paper (Sanchez et. al. in preparation). 

Our ensemble of mock catalogues allowed us to show 
that it is possible to construct an accurate model of the 
shape and amplitude of the REFLEX II power spectrum. 
This model includes the effects of the non-linear evolution 
of density fluctuations, redshift-space distortions and halo 
bias, which introduce deviations in the clustering signal with 
respect to the simple predictions of linear perturbation the- 
ory. This will allow us to use the full information contained 
in the REFLEX II power spectrum to obtain constraints on 
cosmological parameters. 

This paper is organized as follows. In Section [2] we de- 
scribe the REFLEX II sample and the survey selection func- 
tion, followed by a brief description of the construction of 
the mock catalogues in Section ^. 31 In Section[3]we describe 
the power spectrum estimator and show the measurements 
of the REFLEX II window function and the covariance ma- 
trix. In Section 14.11 we model the amplitude of the power 
spectrum measured from the mock catalogues. In Section 
14.21 we explore the sensitivity of the REFLEX II sample to 
distortions induced by flux-selection effects. In Section 14.31 
we model the shape of the power spectrum. The model of 
the shape and the amplitude is applied to the REFLEX II 
sample in Section \S\ We end with our conclusions in Section 

El 

Our fiducial cosmological model consist of a flat ACDM 
Universe with a matter energy density parameter of fi m = 
0.25, a dark energy equation of state w = —1, a dimension- 
less Hubble parameter h — 0.7 and a spectral index of 
primordial scalar fluctuations n s = 1. Throughout this pa- 
per we always refer to the X-ray luminosity in the ROSAT 
hard energy band 0.1 — 2.4 keV and, whenever it is not ex- 
plicitly written, its units are given in 10 44 ergs -1 /i -2 . 



3 The Hubble constant Ho in units of 100 km s 1 Mpc l . 
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Right ascencion a[deg] 

Figure 1. Distribution of REFLEX II clusters in equatorial coordinates. Open circles represent the position of the REFLEX II clusters. 
Filled circles represent the REFLEX II clusters without redshift. The points represent 5 per cent of the random catalogue constructed 
with the REFLEX II selection function. The empty regions in the Southern hemisphere corresponds to the cut in galactic coordinates 
1 6.17 1 < 20° of the Milky Way (band) and the Magellanic clouds. 



2 THE REFLEX II GALAXY CLUSTER 
CATALOGUE 

2.1 Catalogue construction 

Fig. [TJshows the angular positions of the 911 galaxy clusters 
of the REFLEX II sample. The catalogue covers 13924 deg 2 
(4.24 sr) in the southern hemisphere (S < 2.5°) where 
the Milky Way and the Magellanic Clouds have been ex- 
cised to avoid contamination due to stars and regions with 
high X-ray absorbing neutral hydrogen column density and 
high extinction in the optical band. The limiting flux of 
the REFLEX II sample is 1.8 x 10" 12 erg s _1 cm -2 and 
includes 464 new clusters in addition to the original RE- 
FLEX sample, which had 447 clusters to a limiting flux of 

3 x 1CF 12 erg s _1 cm -2 . The detection tech nique is the same 
as that developed for th e REFLEX sample (fBohringc r et alj 
|2001| ;|G uzzo et al]|2009l ) , with spectroscopic redshifts for 860 
clusters and spanning luminosities in the range 4.9 x 10 40 < 
L x /(erg s _1 /i _a ) < 1.96 x 10 45 . The missing « 6 per cent 
of the redshifts will be obtained in observations made at La 
Silla in fall of 2010 and spring 2011. The angular distribu- 
tion of X-ray clusters without redshift in the REFLEX II 
sample is also shown in Fig[TJ (filled circles), displaying no 
particular pattern on the sky. We estimate that our results 
will not change substantially due to this incompleteness. 

While a more detailed description of the sample con- 
struction and the derivation of the cluster parameters will 
be given in a forthcoming paper (Bohringer et al. in prepa- 
ration), we provide here a brief description of these mea- 
surements and calculations. Source counts for the galaxy 
clusters in the RASS have been determined in the 0.5 to 
2 keV energy band by means of the g r owth curve analy- 
sis method described in lBohringer et all l|2000l ). The growth 
curve method is tailored to maximize the aperture in which 
the source counts are determined. The count rate (obtained 
by reference to the exposure maps of the RASS) is then con- 



verted to a nominal flux, F n , by means of XSPEC l|Arnaudl 
Il996h assuming a MEKAL plasma model for a cluster tem- 
perature of 5 keV, a redshift of z = 0, a metallicity of 0.3 
solar and an value for the hydrogen column density taken 
from the measurements of iDickev fc Lockmanl <|l993D . The 
flux limit is imposed on the cluster sample using this value of 
F n . For clusters whose redshifts are known, an improved flux 
value, Fx, is determined by recalculating the count rate to 
flux conversion for a temperature estimat ed via the X-ray lu - 
minosity temperature relation as given bv lPratt et al. I (I2009D 
and by including the proper band corrections (analogous to 
the optical K-correction) for the actual cluster redshift. The 
measured flux is converted to an estimated flux within an 
aperture radius of rsoo (the radius in which the mean matter 
density of the cluster is 500 times the cr itical density of th e 
Universe) by means of relations given in iPratt et all l|2009l) . 
The flux extrapolation (and in some cases interpolation) is 
achieved b y assuming a cluster surface bright ness following 
a /3- model i|Cavaliere fc Fusco-Formiano|[i976r ) with P = 2/3 
and core radius r c = (1/7) rsoo. The cluster rest frame X-ray 
luminosity is calculated by means of the luminosity distance 
from Fx by taking the proper band corrections for the red- 
shift into account. 

For the determination of the sky-position dependent se- 
lection function, the minimum luminosity as a function of 
redshift and the position on the sky L^- m (a, 8, z) is calculated 
assuming the nominal flux limit and taking into account that 
for 5.4 per cent of the sky the exposure is too short (< 100 
sec. mostly due to instrument shut-down during passages 
of the radiation belts in South Atlantic Anomaly) to reach 
the nominal flux limit. The values of Lx m («,<5, z) are then 
derived by accounting for the proper Fx for given redshift 
and by performing an iterative backward engineering of the 
above described process. The surveyed area defining the RE- 
FLEX II sample has been divided into A r p i x = 13902 pixels 
with area ~ 1 deg 2 . For each pixel centered on equatorial co- 
ordinates (cti,5i) the limiting luminosity L^ m («i, Si, z) was 
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Sample i 




Ni 




N V LS 


"VLS 


1 


0.015 


661 


0.070 


32 


4.318 


2 


0.049 


625 


0.080 


90 


1.856 


3 


0.154 


512 


0.099 


150 


0.6562 


1 


0.245 


441 


0.112 


159 


0.3905 


5 


0.490 


306 


0.136 


154 


0.1535 


6 


0.588 


260 


0.143 


157 


0.1155 



Table 1. Minimum X-ray luminosities L™ m (i = 1, • • • , 6) used 
to define sub-samples of the REFLEX II catalogue, with a max- 
imum redshift z max = 0.22. Ni denotes the number of clus- 
ters in each sub-sample and Zi is the corresponding mean red- 
shift. Ny^c; is the number of REFLEX II clusters and fiy^g the 
mean density in each sample volume-limited sample, in units of 
10 _6 (Mpc h -1 ) -3 . X-ray luminosity in units 10 44 erg s -1 h~ 2 in 



the energy band 0.1 — 2.4 keV. 
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Figure 2. Luminosity-redshift diagram for the REFLEX II sam- 
ple (points). Dashed lines schematically represents the volume- 
limited samples defined by the limiting luminosities of Table [JJ 



tabulated in the range ^ z ^ 0.8. Given the minimum 
count rate of 20 cts a nd the geometrica l boun daries of the 
survey (see table 1 of iBohringer et all (1200 ll )) we end up 
with a total of 787 galaxy clusters. 



2.2 The REFLEX II Selection Function 

We measured the REFLEX II X-ray luminosity function 
using the V max estimator. We fitted the measurements by 
means of what we call an extended Schechter function, 



$(L x )dL } 



no 



Lx_ 

Li 



e q (x) 



(1) 



wher e e„(x) deno tes the so-called q— exponential function 
(e.g. lTsalWl2009h defined as 



e q (x) = 



(1 + x(l - q)) 1 ^ 1 -^ 



g = i, 



(2) 



This parameterization provides a better description of 
the high luminosity tail of the X-ray luminosity function 
(Sanchez et al., in prep). We implemented a Markov Chain 
Monte Carlo (MCMC) technique to determine the best fit 
parameters (no, a, L+, q) given by 



no 
a 
Li 



(4.081^:^) 



x 10~ 6 (Mpc/T 



-1.536 



-0.006 
-0.0135 



(0.6283lg"Z) 



x 10 44 erg s 



1 h 2 , and 



= 1.313Z 



(3) 



In order to explore the behavior of the clustering strength 
as a function of X-ray luminosity, we split the REFLEX II 
sample in six sub-samples characterized by a minimum lu- 
minosity L™ in . These are listed in Table [1] together their 
corresponding number of clusters Ni and mean redshift Zi . 
For a given luminosity cut L™ m , the expected number den- 
sity at a position r, characterized by angular coordinates 
(a, S) and redshift z, is obtained via integration of the X- 
ray luminosity function as 



poo 

n(r;Lf n )= / $(L X ) dL 

JL(r) 



(4) 



where the lower integration limit L(r) is given by the RE- 
FLEX II sensitivity map as 



£(r) = Max(Lr ln ,Li m (a,5,«)). 



(5) 



We also constructed six volume-limited samples (hereafter 
VLS) using the same minimum luminosities. These VLS are 
schematically represented in Fig. [2] Table [T] also lists the 
number of clusters, the mean redshift and the cluster number 
density for these sub-samples. 

Finally, we created a random catalogue of N T = 2 X 
10 6 objects with luminosities greater than Lx = 1.4 x 
10 43 erg s _1 h~ 2 . In Fig. \T\ we show the angular positions 
of a subset of the random sample. Note the variations in the 
angular distribution of the random catalogue, which follow 
the fluctuations in the sensitivity map of the REFLEX II 
survey. 



2.3 The mock catalogues 

In this section we give a short outline of the construction 
of a suit of REFLEX II mock catalogues. A more detailed 
description and analysis will be given in a forthcoming paper 
(Sanchez et al. in prep). 

We used the z = output of the 50 realizations from 
the L-BASICC II TV -body simulations (|Angulo et al.ll2008l ; 
ISanchez et~aiT 2008a) to generate a suit of REFLEX II mock 
catalogues. The cosmological model adopted in the simula- 
tions consists of a flat ACDM universe with a matter energy- 
density parameter Sl m = 0.237, a baryon energy-density pa- 
rameter of Sib = 0.046, a dimensionless Hubble parameter 
h = 0.73, a dark energy equation of state w — — 1, and 
an initial matter power spectrum characterized by a scalar 
spectral index n s = 0.954 and normalized to as = 0.77, 
which is in close agreement with the latest constraints on 
cosm ological parameters from CMB and LSS measurements, 
(e.g. ISanchez et al.ll2006|; |Spergel et al.l 2007; Sanchez et al.l 
|2009| - iKomatsu et al.ll201ol v" The low value of as is particu- 
larly important in order to obtain cluster number counts in 
accordance with observations. 

Each of the L-BASICC II simulations follows the dark 
matter distribution using 448 3 particles over a comoving box 



© 2010 RAS, MNRAS 000, [TV?? 



The REFLEX II power spectrum 5 



of side 1.34Gpcft _ . We used halo catalogues identified via 
a friend-of-fricnd (FoF) algorithm with a linking length of 
b — 0.2. The resulting halo mass resolution is M = 1.75 x 
10 13 Mq/T 1 . 

We assigned luminosities to the dark matter haloes in 
the simulations using a mass-luminosity M — Lx relation 
with an intrinsic scatter a\ n L = 0.26. This value was mea- 
sured bv lStanek et al.l (|2006t ) under the assumption of a flat 
cosmological model with f2 m = 0.24, close to our fiducial 
cosmology. Flux errors are included in the luminosity as- 
signment by assuming a fixed error in Lx of SLx / Lx = 20 
per cent, which characterizes the flux errors in the REFLEX 
II sample fairly well. For the mean M — Lx relation we as- 
sumed a power law with a mass-dependent slope. Setting 

1 = l0g10 ( lO^ergs-ife-0 ' m = l0Sl0 ( lO"M fc-O ' 

(6) 

our mass-luminosity relation reads as 

£ = a + bm + cm 2 . (7) 

Using a MCMC technique, we calibrated the set of coeffi- 
cients (a, b, c) such that the resulting luminosity distribution 
of the illuminated dark matter haloes match the measured 
X-ray luminosity function from the REFLEX II sample. The 
best fit values are a = -1.3164, b = 1.8769 and c = -0.2955. 
Changes below 1 per cent in these parameters are observed 
when the flux error is reduced to 10 per cent, although this 
value underestimates the observed flux errors of the RE- 
FLEX II sample. The goodness of this self-calibration is 
confirmed by comparing the REFLEX II luminosity func- 
tion and the mean luminosity function determined from the 
mock catalogues, obtaining differences of ^ 3 per cent. In- 
terestingly, this calibration of the M — Lx relation shows 
a deviation from the power-law behavior at low luminosi- 
ties. This is in agreeme nt with recent observations and N- 
body hydro-simulations (|Puchwein et al . 2008; S tanek et al.l 
2010). Nevertheless we do not attempt to extract conclu- 
sions on the underlying mass-luminosity relation for a num- 
ber of reasons. First, the FoF halo-finder algorithms over- 
estimates the halo mass due to systematics introduced by 
the fi nite number of part icles used to define a dark matter 
halo. IWarren et ail (|200fJ ) showed a way to correct for this 
systematics, which only depends on the FoF mass. For the 
lowest luminosity cut shown in Table [TJ the distribution of 
minimum halo masses selected by equation (0 in the mocks 
peaks at M ~ 5 x 10 1 3 Mp, h -1 , which, following the correc- 
tion of lWarren et al.l |2006f ) corresponds to an offset of ~ 40 
per cent with respect to unbiased mass estimations. Never- 
theless, we do not attempt to correct the L-BASICC II FoF 
masses, for any correction would lead to a new set of pa- 
rameters (a, 6, c) which will still reproduce, by construction, 
the observed X-ray luminosity function. Secondly, there is 
not a clear one-to-one relation bet ween FoF masses and the 
spherical overdensity masses (e.g. iLukic et al.ll201Cl ), which 
are usually im plemented in the calibration of cluster scaling 
relations (e.g. iPratt et al.ll2009l ). 

We next transformed the coordinates of the haloes to 
redshift-space using r — > r + v • r/Ho where v is the pe- 
culiar velocity of the center of mass of the haloes. We ne- 
glected spectroscopic redshift errors for being of the order of 




0.03 0.06 0.09 0.12 0.15 



kj [h Mpc- 1 ] 

Figure 3. Example of the REFLEX II shell averaged window 
(matrix) function Wij for different modes ki determined for the 
sub-sample with limiting luminosity L™ m . Each curve is normal- 
ized to f dkjW(ki, kj) = 1. 

a ~ 10~ 3 (jGuzzo et al.ll2009l ). Finally we observed the illu- 
minated haloes through the REFLEX II mask and obtained 
a set of mock catalogues with the abundance and geometry 
in agreement with those of the REFLEX II sample. 

We constructed two sets of mocks. One set consists of 
50 mocks covering the full volume of the REFLEX II sam- 
ple (out to z ~ 0.5). For the second set, we noticed that 
the effective volume of the REFLEX II sample reaches a 
maximum at redshift z « 0.2 (see Section [3 . 2 p . Therefore, 
including clusters with higher redshift will not help to im- 
prove the signal-to-noise ratio of the measurements. For this 
reason we set a maximum redshift of z — 0.22 which allowed 
us to construct 100 independent mock catalogues out of the 
50 L-BASICC II realizations. Unless otherwise stated, we 
use the set of 100 mocks to carry out our statistical analy- 
sis. The set of 50 mocks were used for consistency checks. 

In the analysis of the REFLEX power spectrum, 
ISchuecker et alj (|200ll ) constructed ten mock catalogues by 
assigning X-ray luminosities to the dark matter haloes in a 
s mall set of A-body simulation s using the M — Lx relation 
of IReiprich fc Bohringerl (|l999l ). No intrinsic scatter around 
the mean M — Lx relation, flux errors or missing flux correc- 
tions were introduced in the mocks. Note that this difference 
with respect to our analysis will result in a different predic- 
tion of the amplitude of the power spectrum for a given 
luminosity cut. In Section 15.21 we address the modeling of 
the amplitude of the power spectrum taking these details 
into account. 



3 METHODS 

3.1 The measurement of P(k) 

We have measured the power spectrum of the REFLEX II 
sample containing clusters with luminosities Lx J? In* 1 *, 
which represents 760 objects in the redshift interval cov- 
ered by the REFLEX II survey. We implemented the stan- 
dard minimum variance w eighting power spectrum estima- 
tor of Fel dman et al.|[l993 (hereafter FKP ), which defines 
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a weighted density fluctuation 
F(r) = w(r) (n c (r) - cm r (r)) , 



(8) 



where n c (r) (n r (r)) is the number density of clusters in the 
real (random) catalogue for a given luminosity cut (where 
we dropped the L mm dependence to avoid clutter). The pa- 
rameter a, given by 



J w(r)n c (r) d 3 r 
J w(r)n r (r) d 3 r 



(9) 



forces the fluctuation to have zero mean, J P(r) d 3 r = 
0. The optimal nor malized weights w(r) are given by 
iFeldman etail (| 19941 ) 



■ fi(r)Pe. 



h 



+ n(r)P ci 



■<Tr 



(10) 



where P GSt is an estimate of the amplitude of the power 
spectrum expected to be measured, for which we have cho- 
sen P es t = 2 x W 4 (Mpch~ 1 ) 3 . We used the FFTw al- 
gorithm (|Frigo fe Johnson! 120051 ) embedding the REFLEX 
II volume in a cube divided in iV gr id = 512 3 cells and 
implemented a triangular s haped cloud mass assignment 
l|Hocknev fe Eastwoodlll988l ) correcting afterwards for alias- 
ing effects. The length of the sides of the cube are determined 
by Lbox = 2r(z max ), with z max = 0.22 which corresponds to 
a box size of 1.25 Gpc h" 1 for our fiducial cosmology. The 
fundamental mode is Afc = 27r/L bo x = 0.0049 h Mpc -1 . The 
Nyquist frequency for this box is A?Ny q = l-27h Mpc -1 , and 
we can ignore aliasing effects on wavenumbers smaller than 
k ~ 0.7 ft Mpc - . We subtract the shot noise and average in 
spherical shells to obtain the bin-averaged power spectrum 
P(k). This measurement is the convolution of the underly- 
ing cluster power spectrum with |iy(k)| 2 , the square of the 
Fourier transform of the REFLEX II window function given 
by 



W(k) = / n(r)ty(r)e" 



"dV 



(11) 



We follow the procedure of ICole et al.l (|2005h to con- 
struct the window function in matri x form Wg, by u sing a 
Gauss-Legendre integration scheme (|Press et alj|2003 ). The 
measured power spectrum can be written as a matrix mul- 
tiplication 



Wu 



(12) 



where P(kj) is the underlying power sp ectrum and the term 
Wip accounts for the integral constraint l|Percival et al.ll2007l ; 
iReid et al] |2010| ). As an example, Fig. [3] shows some ele- 
ments of the window matrix of the sample defined by the 
minimum luminosity L™ 111 (panel a) and the volume lim- 
ited sample defined by the same luminosity cut (panel b). 
As expected, large scale modes receive contributions from 
intermediate and even small (k 0.3h Mpc -1 ) scales. 

We used the window matrix to asses the possibility of 
detecting the signal of the BAO in the measured REFLEX 
II power spectrum. The signature of BAO in the dark mat- 
ter halo distribution of the L-BASICC II simulations has 
been analyzed both in the s patial two-point correlation func- 
tion JSanchez et al l l2008al ) as well as in the power spec- 
trum ^ngul^^y|20ol)^ We used the fitting formulae of 




0.01 0.05 0.1 

Wavenumber k [h Mpc -1 ] 

Figure 4. a) Mean power spectrum from the mock catalogues 
(open points) with its corresponding standard deviation (dotted 
lines), for the sub-sample characterized with the limiting lumi- 
nosity L™ ln . The solid line P n -w(k) represents a non-linear matter 
power spectrum P(k) without BAO and convolved with the RE- 
FLEX II window function. The dashed line represents the same 
theoretical prediction with BAO, also convolved with the win- 
dow function, b) Ratio of the mean mock power spectrum (open 
points), and the power spectrum P(k) to the power spectrum 
fnw(fc). Dotted lines denotes the la standard deviation. 



lEisenstein fe Hrj (|l998l ) to compute the matter power spec- 
trum for the cosmological model of the L-BASICC II sim- 
ulations, i ncluding a non-lin ear correction computed with 
HALOFIT (ISmith et alj 120031 ). We also computed a model 
with the same broad-band shape but with out any baryonic 
oscilla tions, P nw (fc), provided as well by lEisenstein fe Hrj 
l|l998l) . Panel a) of Fig. H shows the comparison of these 
theoretical models (solid and dashed lines), convolved with 
the REFLEX II window function, with the mean power spec- 
trum from the mock catalogues (open circles). It can be seen 
that the difference between these two models is much smaller 
than the standard deviation in P(k) that corresponds to the 
REFLEX II volume, which can be determined from the en- 
semble of mock catalogues (dotted lines, see Section |3~2)) . 
This can be more clearly seen in panel b) of the same fig- 
ure, which shows the ratio of these power spectra to P n w(fc). 
The convolution with the window function damps the acous- 
tic oscillations in the power spectrum to a level where they 
cannot be distinguished from a model without BAO. Fur- 
thermore, we computed the x 2 of these two models (ana- 
lytically marginalizing over the amplitude as described in 
iLewis fe Bridle 1 120021 ) and found a difference of less than 
~ 3 per cent between them. We thus conclude that, due to 
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the survey volume, no statistically significant signal of BAO 
can be detected in the REFLEX II power spectrum. 

A detection of BAO signatures in a galaxy cluster sam- 
ple was claimed by Miller et al. (2001) based on the mea- 
sureme nts of the power sp ectrum of the Abell/ACO cluster 
survey (jMiller et all200ll ). This claim was based on different 
statistical methods from the ones applied in this work and 
no mock catalogues were used to asses the statistical signifi- 
cance of the detection. Our cluster sample contains approxi- 
mately the same number of objects as that used by Mi ller et 
al. (20 01), while probing a larger volume. Recently, iHiitsil 
(2010) reported a 2a detection of BAO in the maxBCG 
galaxy cluster survey, which probes a larger effective vol- 
ume than the REFLEX II sample. 

Finally, regarding the redshift incompleteness of the 
REFLEX II catalogue mentioned in Section 2.1, we veri- 
fied that our results are not substantially modified when 
the power spectrum is measured after randomly subtracting 
up to 20 per cent of the total number of clusters. 

3.2 Covariance matrix 

The determination of the covariance matrix of the power 
spectrum is a key step towards obtaining constraints on cos- 
mological parameters. Its determination represents a non- 
trivial task both from theoretical and numerical perspec- 
tives. From an analytical point of view, the covariance ma- 
trix can be decomposed into a Gaussian component, which 
depends on the measured power spectrum, and a non- 
Gaussian contribution related to the bi- and tri-spectrum 
ijMatarrese et al.ll 19971 : IVerde fc Heavendl200ll : ISmithll2008h . 
On large scales, which are well described by linear pertur- 
bation theory, these non-Gaussian contributions vanish. On 
small scales however the non-linear evolution introduces cou- 
pling between Fourier modes generating signatures of non- 
Gaussianities. Hence, to precisely determine the covariance 
matrix, it is necessary to have a model for these high-order 
statistics, together with redshift distorti ons, survey window 
function effects, beat-coup ling effects l)Rimes fc Hamilton! 
120061 ; ITakavashi et al. 20091') and corr elations introduced by 
bin-averages jMeiksin fc White! 1 19991) . From the numerical 
side, recent experime nts with iV— body simulations (e.g. 
ITakavashi et al.l 120091 ) have shown that large numbers of 
mock catalogues (or realizations) are required for a precise 
determination of the covariance matrix as needed to con- 
strain cosmological parameters. 

We used our ensemble of 100 mock catalogues to obtain 
an estimate of the bin-averaged covariance matrix C(ki, kj) 
of the REFLEX II power spectrum by 



AT, 



pn 



(13) 



where P™ = P n (ki) is the measured power spectrum in the 
n-th mock catalogue in the bin centered at ki and Pi is the 
mean power spectrum from the ensemble of mocks at the 
same bin. As an example, panel a) of Fig. shows examples 
of the correlation coefficients defined from the covariance 
matrix via f%j = Cij /(CaCjj) 1 ^ 2 for L™ 1 " (upper triangu- 
lar part) and I/™ 111 (lower triangular part). For comparison, 
panel c shows the correlation matrix inferred from the clus- 
tering of the illuminated halos in the L-BASICC II simu- 
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Figure 6. Effective volume probed by the REFLEX II sample 
(see equation 1156 as a function of the maximum redshift for three 
different cuts in luminosity. 



lation for the same luminosity cuts. The covariance matrix 
of the mocks contains important off-diagonal terms which 
arise from the mode coupling induced by the window func- 
tion. All the statistical analyzes performed in this work are 
bas ed on the covariance m atrix defined by equation (|13[) . 

iFeldman et al.l (Il994f ) derived an approximated expres- 
sion for the variance of the spherically averaged power spec- 
trum under the assumption that the Fourier modes are 
Gaussian-distributed. This is given by 



a 2 {k) _ 2 
P(fc) 2 ~ V k V eff (k)' 



(14) 



where Vk ~ 47rfc 2 <5fc/(27r) 3 is the volume of a spherical shell 
of width 8k and V c g(k) is the effective (co herence) volume 
probed by the survey at a scale k, defined bv lTegmarkl (1997) 
as 



V eS (k) 



n(r)P(k) 
1 + P(k)n(r) 



(15) 



Equation (|14[) assumes that the power spectrum P(k) is 
smooth on scales 5k and applies for wavenumbers fc< S> Sk. 
Fig. [6] shows the effective volume for different luminosity 
cuts as a function of the maximum redshift of the sample. 
As pointed out in Section 12.31 the effective volume gained 
by including objects with redshifts z 0.22 is very small 
and only leads to an increase of the shot-noise. 

Fig. [7] shows a comparison of the theoretical variance 
computed using equation (|14fl with that derived from the di- 
agonal elements of the covariance matrix of the mocks (equa- 
tion Q3J for the luminosity cuts defined in Section [2.21 The 
theoretical predictions were computed using a linear theory 
power spectrum with an amplitude rescaled to match that of 
the REFLEX II measurements for the correspondent lumi- 
nosity cut. This comparison shows a very good agreement 
between the results obtained from the ensemble of mocks 
and the theoretical prediction for all values of L mm . Equa- 
tion (|14|l thus provides a good estimate of the error bars in 
the measured power spectra. 
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Figure 5. Comparison of the correlation matrix \r(ki,kj)\ of the power spectra measured with a) the FKP estimator and b) the PVP 
estimator (see Section 14. 2H for two luminosity-cuts, c) Correlation matrix determined from the illuminated dark matter halos of the 
L-BASICC II simulation with same luminosity-cuts. 
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Figure 7. Bin-averaged standard deviation a p (k): solid line rep- 
resents the standard deviatation determined from the set of mock 
catalogues following equation (1 1 3 D - The dotted line shows the pre- 
diction form the equation 1141 . 



4 UNDERSTANDING THE REFLEX II 
POWER SPECTRUM 

In this section we use our ensemble of mock catalogues to 
understand what can be expected from a measurement of 



the power spectrum in the REFLEX II catalogue. Section 
4. II deals with a theoretical model for the luminosity depen- 
dence of the bias. In Section 14.21 we analyze the possible 
presence of a scale-dependent systematic bias induced by 
the flux-limited selection. Finally, Section 14.31 describes a 
model for the full shape of the power spectrum including 
the distortions introduced by non-linear evolution. 



4.1 A theoretical prediction of luminosity bias 

Assuming that on the scales of interest the underlying 
halo-mass bias is scale-independent, the power spectrum of 
galaxy clusters of a given luminosity Lx and at a given red- 
shift 2 can be written as 



P c i(fc, z; Lx) = b (L x ,z)P nmt (k,z), 



(16) 



where P ma t (fe, z) is the matter power spectrum (we now drop 
the redshift dependence of these factors for simplicity). The 
luminosity bias b(Lx) is given in terms of the halo mass 
func tion n(M) an d the underlying halo-mass bias b(M) as 
(e.g. ICooravll2006l ) 



b{L x ) = 



/ n{M)b{M)p{L x \M) dM 
/ n(M)p(L x \M)dM 



(17) 



Here p(Lx\M) represents the conditional probability distri- 
bution of assigning a X-ray luminosity Lx to a dark matter 
halo of mass M. As described in Section 12.31 for our mock 
catalogues p(Lx\M) is given by a log-normal distribution 
with a mean given by equation © and a scatter including 
an intrinsic dispersion <7i n L = 0.26 and a 20 per cent lu- 
minosity error added in quadrature. In equation (|17p . the 
numerator is the theoretical prediction for the X-ray lumi- 
nosity function 



9(L, 



n(M)p(L x \M) AM, 



(18) 



given by the halo mass function and the mass-luminosity 
relation. 

In order to test the predictions of equation iJTTJ, we 
used the 50 realizations of the L-BASICC II to measure the 
absolute luminosity bias. We assigned luminosities to the 
dark matter haloes using the calibrated mass-luminosity re- 
lation described in Section \2. 31 We then measured the halo 
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Figure 8. Measured luminosity bias b(ALx) from the L-BASICC 
II simulation in real space (filled points) and redshift space (open 
points). Solid and dashed lines represent the prediction from equa- 
tions I I17II and ( 12 1 1) respectively. 



power spectrum P c \(k; ALx), in real-space, by splitting the 
sample in luminosity bins of width ALx ■ We used the mea- 
surements of the real-space dark matter power spectrum of 
the same simulations to determine the ratios 



b(ALx 



P c i(k,AL> 



(19) 



Pmat(fc) 

For this we used Fourier modes in the range 0.01 ^ 
k/(h Mpc -1 ) ^ 0.1, where this ratio is compatible with a 
scale-independent bias factor. The obtained measurements 
are shown by the filled circles in Fig. [8] 

The measured bias factors can be estimated by inte- 
grating equation (|17[) within the given luminosity bin 



b(ALj 



f ALx $(Lx)b(Lx)dL x 



(20) 



where <fr(Lx) is given by equation (|18|l . This prediction is 
shown by the solid line in Fig. [5] F or this we used the mass 
function from Jenkins et al.l (|200ll ) and the halo bias from 
iTinker et al.l (|2008l ), which provide an excellent description 
of these quantities in the L-BASICC II simulations. There 
is an excellent agreement between the theoretical b(ALx) 
and the measurements obtained from the simulations. This 
agreement depends on the accuracy with which the M — Lx 
relation is known. 

Under the assumptions of linear evolution and the dis- 
tant observer approximation, the cluster power spectrum in 
redshift-space is a boosted version of its real -space counter - 
part, P^(k;AL x ) = S(ALx)Pci(k;AL x ) dKaiser 1 1 19871) . 
where 



S(AL> 



1 + 



/ 



b{AL x ] 
is the f 



(21) 

rowth index (e.g. 



3 6(AL X ) 5 
where / = d In D (a) /d in a 

■Wang fc SteinhardtJ Q9981 ) and D(a) is the growth factor 
(e.g. .Peebles! Il980l ). The open circles in Fig. [8] show the 
bias factors 6 s (ALx) estimated from the L-BASICC II sim- 
ulations as in equation (|19[) using the redshift-space halo 
power spectrum P^(k, ALx)- The dashed line in this figure 
corresponds to the theoretical prediction for this quantity, 
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Figure 9. Effective luminosity bias measured from the REFLEX 
II mock catalogues (filled circles with error bars). The solid line 
shows the prediction from equation I I25H with the fiducial values 
a lnL = 0.26 and constant luminosity error of 20 per cent. The 
dashed line is the prediction for a VLS. 



given by b(ALx) 2 S(ALx)- The agreement between this pre- 
diction and the results from the A^-body simulations show 
the validity of this simple treatment of the redshift space- 
distortions. 

The results of equations (|20|l and (|21[) can also be used 
to estimate the effective bias of a power spectrum measure- 
ment in the REFLEX II catalogue for a given luminosity 
cut L mm . In order to achieve this, we first write the clus- 
ter power spectrum in redshift-space at a given redshift z in 
terms of the linear dark-matter power spectrum at z = as 



P cl (fc, z\ > L 
where 



x\ 7 s / r mm \ 2 

) = 6 (> L ,z) 



t(M = 0), 



b E {z, > L mm ) 2 = b(> L mm , z) 2 S{z, > L mln )D 2 (z). 



(22) 



(23) 



In this expression we take into account that the minimum 
luminosity included in the sample varies with z following the 
REFLEX II selection function. Accordingly, we determine 
the bias b(z, > L mm ) following equation (|20[) with the red- 
shift dependence given by the REFLEX II sensitivity map: 



b(z, > L min ) = - 



n(z, > L min ) J L(Z) 



$(L)b(L) dL, 



(24) 



where the mean number density n(z, > L mln ) and the lower 
integration limit L(z) are given by equations @ and (0 
respectively. Note that L(z) should not only depend on the 
redshift but also on the angular position, according to the 
REFLEX II sensitivity map. However, in order to give an 
estimate of the effective bias at a fixed redshift we made 
an average of the values of L mm (r) of all the A^pix pixels 
within the REFLEX II mask. Individual pixels displayed 
small differences compared to the average in the final result. 

The effective bias of the full sample will then be given 
by the average of the bi as factors of equation l|24p over the 
observed volume as (e.g. ISuto et al]|2000l ; iMoscardini et al.l 
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Figure 10. Mean squared luminosity bias for pairs separated by 
a scale r in the REFLEX II mock catalogues, for the sub-sample 
L™ m (filled circles) and its corresponding VLS (open circles). 



120001 ) 



feeff(> L ml T = 



„iz, > L min ) 2 ^dz 



(25) 



where the integrals are evaluated in the redshift interval 
of the sample with a volume element dV/dz = r(z) 2 / H(z) 
according to our fiducial cosmology. This therefore allows us 
to write a prediction for the observed power spectrum from 
a sub-sample characterized by the luminosity cut L mm as 

P c i(fc; > L min ) = 6 cff (> L min ) 2 P mat (fc, z = 0). (26) 

The solid line in Fig. [9] shows the prediction for the effec- 
tive bias of the REFLEX II catalogue as a function of the 
minimum luminosity L mm computed using equation (|25|) 
(setting D(z) = 1 since, by construction, the mock cata- 
logues assume no redshift evolution). For comparison, the 
dashed line in Fig. [9] represents the equivalent prediction 
for a volume limited sample. Equation ()25|) gives an excel- 
lent description of the direct measurements obtained from 
the mock catalogues (shown by the filled circles). Then, this 
model provides a means to extract the important cosmolog- 
ical information contained in the amplitude of the measured 
REFLEX II power spectrum. 

4.2 Systematics: distortions induced in a flux 
limited sample 

Due to the flux-limited nature of the REFLEX II catalogue, 
large scales might be probed predominantly by the most lu- 
minous clusters, with a higher clustering amplitude. This 
would artificially increase the measured power spectrum on 
large scales, introducing a scale-de pendent distortion with 
respect to the volume limited case. ISchuecker et alj (|200ll ) 
analyzed this problem in detail for the REFLEX sample 
and concluded that no significant effect can be detected for 
scales r < 150 Mpc h . In this section we perform a sim- 
ilar analysis on the REFLEX II sample. As the REFLEX 
II sample spans a wider range of luminosities, an d cove rs a 
larger volume than that used by ISchuecker et al] (|200ll ). it 
is necessary to test whether this systematic effect can affect 
our measurements. 



Fig. [10] shows the number of pairs with separation r 
weighted by the individual biasing factors of the pair mem- 
bers (computed using equation l|17[)) determined from the 
mock catalogues, yielding the average squared bias factors 

(* ! W) = iE^)Hii)- (27) 



n(r) 



where the sum is done over pairs separated by scales in 



the range 



< \n — rj 



< r + iA and n(r) is the 



number of clusters in the same interval. The closed circles 
show the results obtained from the sample with minimum 
luminosity L™ m (see Table [TJ , while the open circles cor- 
respond to the equivalent measurement from the VLS de- 
termined by the same minimum luminosity. The error bars 
are drawn from the standard deviation of the ensemble. 
Two prominent features can be observed from Fig. 1101 On 
one hand, there is an increase in the mean bias on scales 
r w 10Mpc/i -1 . This is understood as to show that pairs 
of clusters separated by these scales have at least one clus- 
ter with luminosity bias higher than the mean value of the 
sample. This is therefore a consequence of gravitational clus- 
tering. On scales smaller than r ~ QMpc/i" 1 the mean 
bias decreases as a consequence of the halo-exclusion (e.g. 
iPorciani et all 1 199a ). This feature is more evident in the 
sub-sample L™ ln than in its corresponding volume-limited 
counterpart, due basically to the low volume sampled by 
the latter. On the other hand, it can be clearly seen that, on 
scales r > 150 Mpc ft -1 , a systematic increase of the average 
squared bias factor exists in the flux-limited sample com- 
pared with the volume-limited case. This is a direct conse- 
quence of the flux-limited nature of the survey. Naively, the 
scale where this flux-selection distortions is relevant would 
correspond to k = 27t/r = 0.04 h Mpc -1 , suggesting that for 
wavenumbers larger than this limit no scale-dependent dis- 
tortion affects the measurements from the REFLEX II cat- 
alogue. Comparing with other sub-samples, we observe that 
this distortion is damped in the high luminosity cuts, which 
almost behave like volume- limited samples (see Fig. [2}. 

To analyze this issue in more detail we used the FKP 
estimator as implemented bv lPercival et all 2004 (hereafter 
PVP). This is also a minimal variance weighting estimator 
which takes into account the absolute luminosity bias (or 
the relative bias to some luminosity L) to obtain an estimate 
of the power spectrum of the underlying matter (or corre- 
spondingly, the power spectrum of objects with luminosity 
L), free of the distortions ind uced by th e flux-l imited selec- 
tion of the sample. Recently, ICai et all l|2010l ) generalized 
this estimator by introducing a mass-dependent weighting 
scheme that minimizes the stochasticity between the clus- 
ter (or galaxy) field and the underlying matter distribution. 
For the goal of this section it is enough to use the PVP 
estimator. 

Panel a) in Fig. [TTJ shows a comparison of the mean 
power spectra obtained by the FKP (dot-dashed lines) and 
PVP (solid lines) algorithms for the different luminosity cuts 
defined in Table [TJ in our ensemble of mock catalogues. In 
the PVP method, we weighted each object by the inverse 
of its luminosity bias, computed using a fit to the redshift- 
space results shown in Fig. [8] The shaded area in Fig. [TT] 
represents the range of scales used to measure these bias 
factors (0.01 ^ fc(/iMpc _1 ) ^ 0.1). This produces a power 
spectrum normalized as that of the dark matter distribu- 
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Figure 11. a) Comparison of the resulting mean power spectrum 
from the mock sample using the PVP estimator (solid lines) and 
the FKP result (dot-dashed lines) for our six cuts in luminosity 
defined in Table[T] The dashed line represents the measurement of 
the dark matter power spectrum of the L-BASICC II simulations. 
The shaded region shows the scales where the absolute bias is 
fitted (0.01 ^ k (/iMpc -1 ) ^ 0.1) b) Ratio of FKP power spectra 
to the PVP measurements for the six luminosity cuts. 



tion (shown by the dashed line in Fig. Hip . Panel b) shows 
the ratios between the mean power spectra obtained using 
the FKP and PVP estimators for the each luminosity cut. 
In the absence of a scale-dependent distortion, these ratios 
should correspond to the bias factors b^ B (> £min), shown 
by the dotted lines. These ratios show no clear signature 
of a scale-dependent distortion for k) ^ 0.02 ftMpc -1 . Con- 
trary to what might be expected from this systematic effect, 
Fig. [TT] shows a weak indication of a decrease in power on 
large scales for the lower luminosity cuts, which is smaller 
than the standard deviation of the measurements. From this 
analysis we conclude that no significant distortion is intro- 
duced in the shape of the power spectra estimated with the 
FKP method. Note however that this statement applies ex- 
clusively to the REFLEX II catalogu e, as such an effect has 
been detected in ga laxy surveys (e.g. iTegmark et al. 1120061 ; 
IPercival et al.ll2~007l ') . Given the lack of systematic distortions 
in the FKP measurements, we chose to use the FKP esti- 
mator to analyze the data from the REFLEX II catalogue. 
Panel b) of Fig.[5]shows the correlation matrix inferred from 
the ensemble of mock catalogues of the power spectra for 
two luminosity cuts obtained using the PVP estimator. A 
comparison with panel a) shows that the PVP estimator 
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Figure 12. Q-model description of the mock power spectrum. 
The shaded region shows the scale where the Q-model was used 
to fit the mean power spectrum from the mock catalogues. The 
bottom panel show the ratio of the difference between the Q- 
model and the measurements to the standard deviation from the 
mocks. 



induces slightly higher correlations between Fourier modes 
compared to the FKP case. 



4.3 Modeling the shape of P(k) 

In this section we use our ensemble of mock catalogues to 
test a model of the shape of the REFLEX II power spectrum. 
We focus on the clusters with luminosities greater than L™ m . 
The filled circles in Fig. [12] show the mean redshift-space 
power spectrum of the mocks for this luminosity cut with 
error-bars determined from the standard deviation of the 
ensemble. It can be clearly seen that this measurement ex- 
hibits an excess of power at small scales with respect to the 
predictions from linear perturbation theory, shown by the 
dashed line. This is due to the combined effect of non-linear 
evolution and redshift-space distortions. 

In recent years the distortions in the shape of the 
power spectrum produced by these effects have been in- 
tensively studied using large A^-body simulat ions and re- 
cent a dvances in perturbati o n the o ry (e.g. ISmith et al.l 



20071; ICrocce fc Scoccimarrol 120061; lAngulo et all I200I 



Sanchez et alj 2*008al ; iMontesano et all 2010h . These analy- 



ses have produced accurate descriptions of these distortions 
to the level of accuracy demanded by forthcoming surveys, 
which will probe volumes much larger than present day cata- 
logues. Due to the moderate volume probed by the REFLEX 
II catalogue (in comparison to the volume probed by cur- 
rent galaxy redshift surveys), percent-level accuracies in the 

treatment of these effects are not required. 

We now te st whether the Q-mo del of lCole et all (|2005l ) 
(modified as in lSanchez et al. I l2008al ) can provide a good de- 
scription of the non-linearities observed in the power spectra 
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Figure 13. REFLEX II power spectrum (filled circles with error bars) for clusters with luminosities Lx > L mm . The REFLEX power 
spectrum is shown by the open triangles. The error bars for these two measurements are taken from equation JH). For comparison 
we also show the measured power spectrum from the 2df GRS taken from I Cole et al.l J2005l ~) (open circles) . The solid and dashed line 
represent the ACDM power spectrum convolved with the REFLEX II and the 2dfGRS window function respectively, and adjusted to 
match the corresponding spectra. Error-bars exceeding the range of the plot are represented by arrows. 



of our mock catalogues. In this model the shape of the clus- 
ter power spectrum is given by 

1.2 

\2 



P c i(fe, > L) = b e s(> L) 



1 + 

1 + Ak + Bk- 



plin 

*■ mat 



(k), (28) 



where P m at(&) is the linear theory matter power spectrum. 

Although this model was originally developed and cali- 
brated to to describe the power spectrum of the 2dFRGS, its 
applicatio n has been extended t o the analysis of other sam - 
ples (e.g. iTegmark et al. 1120061: iPadmanabhan et aTll2007l ). 
In particular. ISanchez et all l~ i|2008l ) showed that this model 
can give a good description of the clustering of the LRG sam- 
ple from SDSS even though it was not specifically designed 
to do so. At the same time this model does not give a good 
description of the shape of P(k) for the main galaxy sample 
in SDSS. The results from the application of the Q-model to 
A^-body simulations show that it can correctly describe the 
clustering of dark mat ter halos above a given mass threshold 
ITegmark et al. Il200rj ). 



We follow lCole et al l (|2005h and fix the value of A = 1.4 
as obtained from the analysis of A^-body simulations, while 
Q and B are left as free parameters whose values will de- 
pend on the limiting luminosity of the sample. We assumed 
all the cosmological parameters to be known and fitted for 
Q and B ma rginalyzing analyticall y over the amplitude (as 
described in iLewis fc Bridle ! [20021 ) . From this analysis we 
obtain the values Q = 24.9 ± 1.1 and B = 12.0 ± 2.1, cor- 
responding to the sub-sample defined by I/™ 1 "- The best fit 
model obtained this way is shown by the solid line in Fig. 1121 
It can be clearly seen that the model of equation 1)281) gives 
an accurate description of the shape of the mean power spec- 
trum from our ensemble of mock catalogues. This can be also 
seen in panel b) of the same figure, where we show the ratio 
between the difference of the mean mock power spectrum 
and the best fit-model to the variance from the ensemble. 
The parameters B and Q fitting the power spectrum of the 
sub-sample L™ 1 " follow a degeneracy that can be described 
approximately by B(Q) = 0.805Q — 8.15. This degeneracy is 
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Figure 14. Measured REFLEX II power spectrum for the different sub-samples defined in Table [JJ FFTw is carried out in a box of 
side L box = 1263.8MPC/1- 1 and P cst = 2 X 10 4 (Mpc/i- 1 ) 3 . The fundamental mode is 8k = 2?r/L box = 0.0049 h Mpc -1 and the Nyquist 
frequency is kj^ = 1.27?iMpc . Points represent the REFLEX II measurements with error bars drawn from equation II 14 1 . The shaded 
region represents the la standard deviation determined from the mocks catalogues. The solid line represents the mean mock power 
spectrum. 



maintained if the amplitude of the model is fixed according 
to equation (|25[) . We can thus use this degeneracy to reduce 
the number of degrees of freedom when constraining cos- 
mological parameters using the measured power spectrum. 
The best fit value of Q increases with the limiting luminos- 
ity of the sample, varying from Q = 20.7 ± 0.9 for L™ m to 
Q = 44.9 ± 2.3 for L£ in . The general trend in the degen- 
eracy B(Q) is maintained for different luminosity cuts. In 
Section 15.31 we compare the predictions of this model with 
the measurement of the REFLEX II power spectrum. We 
have explictly tested that, by adopting a different value of 
A the best fit value for the parameter Q, and the degener- 



acy between Q and B, are slightly changed, while providing 
equally good fits to the data. 



5 ANALYSIS OF THE REFLEX II POWER 
SPECTRUM 

5.1 Measurements 

The measured power spectrum for the REFLEX II sample 
with limiting luminosity L™'" is shown by the filled points 
in Fig. [13] with error bars drawn from the FKP method 
(see equation [14]) . The solid line represents a ACDM linear 
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power spectrum, convolved with the window function of the 
survey. This theor etical prediction was com puted using the 
fitting formulae of lEisenstein fc HuI i|l998h . with amplitude 
rescaled to match that of the REFLEX II measurement. 
This simple exercise shows that the shape of the REFLEX 
II power spectrum is consistent with the predictions of the 
ACDM cosmological model. 

Fig. [T3] also shows a new estimation of the power spec- 
trum of the original REFLEX sample (open triangles). The 
REFLEX power spectrum has a higher amplitude, as ex- 
pected from the higher flux-limit of this sample, and its 
shape is in good agreement with that of the REFLEX II 
measurement. The larger volume probed by the new cata- 
logue reduces the impact of cosmic variance on large scales, 
where the REFLEX II power spectrum exhibits a higher 
amplitude than the measurement in the original REFLEX 
sample. 

Fig.[l3]also sho ws the galaxy po wer spectrum measured 
from the 2dFGRS (|Cole et al.ll2005h . The dashed line rep- 
resents the same ACDM power spectrum described above 
convolved with the 2dFGRS window function. This shows 
that, once their respective window functions have been taken 
into account, the large-scale (fc < 0.1 ftMpc -1 ) shape of the 
power spectra inferred from the REFLEX II and the 2dF- 
GRS are in good agreement and can be described with the 
same cosmological model. At smaller scales, redshift-space 
distortions and non-linear evolution produce deviations in 
the shapes of these power spectra. 

Fig. [14] shows a comparison of the measured REFLEX 
II power spectra (points with error bars) for the six cuts in 
luminosity described in Section [2.21 and the corresponding 
mean power spectra from the mock catalogues (solid lines), 
with their corresponding la standard deviation (shaded re- 
gions). The error bars of the REFLEX II power spectrum 
correspond to the theoretical prediction of the FKP method 
(see section I3T2")) . We observe that the spectra measured in 
the mocks are compatible within la with the REFLEX II 
clustering up to fc ~ 0.3 h Mpc -1 for all luminosity cuts. 
Notice that the mocks were only calibrated to follow the 
X-ray luminosity function of the REFLEX II sample. This 
agreement allows us to use the covariance matrices inferred 
from the ensemble of mock catalogues when analyzing the 
REFLEX II measurements. 

On small scales the power spectra inferred from the 
mock catalogues are affected by the halo exclusion effect. 
Dark matter halos in the simulations have been counted as 
separate entities when they did not overlap with their radii 
of rFoF- On the other hand, in the REFLEX II catalogue 
clusters have been treated as distinct if their X-ray emission 
does not significantly overlap. Due to the short exposures in 
the RASS, the outer boundary of the X-ray emission (in two 
dimensional images, which is significantly smaller than the 
aperture radius determined from one-dimensional profiles) 
is smaller than the radii of rsoo. This produces differences 
between the REFLEX II power spectra and the results from 
the mock catalogues on scales fc > 0.2 liMpc -1 . 

Regarding the redshift incompleteness of the REFLEX 
II catalogue (around 10 per cent), we verified that our results 
are not substantially modified when the power spectrum is 
measured after randomly subtracting up to 20 per cent of 
the total number of clusters. 
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Figure 15. REFLEX II power spectrum for three different min- 
imum luminosities. Lines represent a linear-perturbation theory 
power spectrum with our fiducial cosmology and amplitude given 
by equation ((25} • Error bars correspond to the standard deviation 
drawn from the mock catalogues. Arrows denote the error bars 
exceeding the range of the plot. 



5.2 Amplitude of the REFLEX II P(k) 

Fig. [15] shows the measurements of the REFLEX II power 
spectra for three of the sub-samples defined in Table [1] The 
increase in the amplitude with increasing minimum luminos- 
ity can be clearly seen, showing the signature of luminosity 
bias. 

In Section T4. II we showed that the measurements of the 
effective bias of the REFLEX II mock catalogues are well de- 
scribed by the predictions of equation (|25[) . In this section 
we confront this prediction with the power spectra measured 
from the REFLEX II sample. In order to avoid using the un- 
derlying dark matter power spectrum, we test equation (|25l) 
by means of the a relative luminosity bias r(Lx), defined as 
the ratio between the power spectrum of the sub-sample de- 
fined by a minimum luminosity Lx to that of clusters with 
luminosities greater than a reference value Lx'- 



r(Lx) = 



b c ff(> Lx) 

&eff(> Lx) 



(29) 



The results are shown in Fig. 1161 The open squares show 
the measurements from the REFLEX II data and the filled 
points correspond to the respective measurement from the 
mock catalogues. The solid line shows the prediction from 
equation (|25[) . while the dashed line is the prediction of the 
effective bias for a VLS. The prediction from equation (25) 
provides a good description of the REFLEX II measure- 
ments in the low luminosity cuts. For the last two luminosity 
cuts the agreement is not as good, although the measured 
bias factors are consistent with the theoretical prediction 
within the errors. As can be seen in Fig. 1141 the spectra 
measured from these luminosity cuts show a low amplitude 
at scales fc ~ 0.1/i/Mpc. This dip in the power spectra, which 
is consistent with cosmic variance, explains why the corre- 
spondent bias measurements in Figure 16 lie below the theo- 
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equation J29D . 



retical prediction. The agreement is improved when the bias 
factors are measured in the range 0.01 < fc(/i/Mpc) < 0.05, 
as is shown by the open triangles in Figure 16. These results 
confirm the validity of equation [25] to model the amplitude 
of the REFLEX II power spectrum. 

Note, finally, that as the theoretical predictions for the 
bias are based on the mass- X-ray luminosity relation cali- 
brated from the mock catalogues, the good agreement found 
with the bias factors inferred from the data suggests that 
for the luminosity range we have explored, the calibrated 
scaling relation provides also a good description of the real 
underlying mass-X ray luminosity relation. 



5.3 Shape of the REFLEX II Pik) 

We used the Q-model to analyze the shape of the REFLEX 
II power spectrum for the sub-sample defined by I/™ 1 ", fol- 
lowing the same procedure as in Section \4. 31 The best-fit of 
the Q-model is shown by the solid line in the upper panel 
of Fig. [17] For this measurement we find Q = 24.7 ± 1.5 
and B = 8.6 ± 1.1, with a degeneracy described by B(Q) — 
0.72Q — 9.25. This degeneracy is maintained when fixing the 
amplitude of the model power spectrum according to equa- 
tion (|25p . As in Fig. 1121 the bottom panel shows the ratio 
of the difference between the model and the measurements 
to the standard deviation determined from the mock cata- 
logues. This shows that a model including a correction for 
non-linearities provided by equation (|28|1 gives a better de- 
scription of the shape of the REFLEX II power spectrum 
than the predictions from linear perturbation theory. 



6 CONCLUSIONS 

In this paper we presented the measurement and analy- 
sis of the power spectrum from the new REFLEX II cat- 
alogue, which is an extension of the original REFLEX sam- 



ple l|Bohringer et alj l200ll ) to a lower limiting flux (1.8 x 
10~ 12 erg s" 1 cm" 2 ). The new sample contains 911 X-ray de- 
tected galaxy clusters of which 860 have measured redshifts 
in the range ^ z < 0.6 and X-ray luminosities in the 
range 4.9 x 10 40 «S L x /(erg s' 1 ^ 2 ) sC 1.96 x 10 45 . The 
total flux and X-ray luminosities are estimated using the 
up-to-date scaling relations based on the REXCESS Survey 
l|Bohringer et al.ll2007l : iPratt et al.ll2009l ). 

The new sample allowed us to perform a detailed study 
of the full shape and amplitude of the power spectrum of 
X-ray detected galaxy clusters. We complemented this anal- 
ysis by using a set of 100 independent mock catalogues con- 
structed to match the selection function of the REFLEX II 
survey. The clustering properties of these mock catalogues 
are in good agreement with those measured in the REFLEX 
II sample. Thus, this ensemble provides a reliable tool to test 
the statistical methods applied to the data. In particular, we 
used the mock catalogues to test a model for the luminosity 
dependence of bias, to construct covariance matrices of the 
the REFLEX II power spectrum and to analyze the possible 
systematic effects that might affect this measurement. 

Due to the flux-limited selection of the REFLEX II sur- 
vey, the clustering pattern of galaxy clusters might be af- 
fected by scale-depend ent distortion, as has been observed 
in ga laxy surveys (e.g. iTegmark et al. 1120061 ; iPercival et alj 
120071 ). Using the mock catalogues, we have shown that these 
distortions might affect the clustering in configuration space 
(i.e., when measured with the cluster-correlation function) 
on scales r ^ 150 Mpc /i -1 , which would naively correspond 
to scales k ^ 0.04/iMpc -1 in Fourier space. In order to test 
the impact of this flux-selection effect on the final measure- 
ments of power spectru m, we implem e nted the luminosity 
dependent estimator of IPercival et al] (|2004h ■ which is de- 
signed to correct for this distortion. We observed that the 
shape of the power spectrum measured by means of the FKP 
estimator does not show significant distortions compared to 
the results from the PVP estimator. This implies that the 
flux-selection of the REFLEX II sample does not introduce 
a significant systematic effect in the measurement of the 
power spectrum of this catalogue. 

The shape of the mean power spectrum from our en- 
semble of mock catalogues is in good agreement with the 
measured power spectrum from the REFLEX II sample, 
and is statistically distinguishable from the linear pertur- 
bation theory predictions on intermediate scales. This im- 
plies a clear signature of non-linear evolution in the X-ray 
cluster spatial distribution. Nevertheless, given the level of 
accuracy of the measurements of power spectrum in the 
REFLEX II sample, it is sufficient to mode l these distor- 
tions using the Q-model of I Cole et al. | (|2005h . We find that 
this prescription provides a good description of the mea- 
surements from the mock catalogues on intermediate scales 
(0.02 ^ k /(/iMpc -1 ) ^ 0.25). This model can also be used 
to describe the shape of the measured REFLEX II power 
spectrum, providing a valuable tool to extract the cosmo- 
logical information contained in the shape of this statis- 
tic. The next generation of X-ray galaxy clusters surveys, 
such as eROSIT^ and WFXtH will provide measurements 



4 http:/ /www. mpe.mpg.de/heg/www/Projects/EROSITA/main. html 

5 http://wfxt.pha.jhu.edu/ 
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Figure 17. Best-fitting Q-model for the REFLEX II mock power 
spectrum (points with error bars). See Fig, H2l for description. 



of the two-point statistics of the cluster population with 
higher accuracy than present-day samples, for which a more 
detailed modelling of non-l i nearities will be requir ed (e.g. 
ICrocce fe Scoccimarrol |2006| ; iMontesano et alj|2010h . 

Our measurements of the REFLEX II power spectrum 
are compatible with the prediction of the ACDM cosmo- 
logical model and shows good agree ment with the previou s 
results from the REFLEX sample (|Schuecker et al.ll200ll h 
save the expected differences due to the lower limiting flux 
of the REFLEX II sample. We showed that our measure- 
ments cannot provide a statistically significant detection of 
BAO, which is mainly due to the moderate volume probed 
by the survey (compared to the volume probed by current 
galaxy redshift surveys). We found that the power spectra 
measured from the REFLEX II sample and the mock cata- 
logues are compatible with a scale-independent effective bias 
in the range of wavenumbers 0.01 ^ k /(hMpc -1 ) ^ 0.1, and 
that a simple theoretical prediction, based on the halo-mass 
bias, the halo mass function and the mass- luminosity rela- 
tion, is able to describe these measurements. This, together 
with the modeling of the shape of the power spectrum given 
by the Q-model, provides a link to the cosmological models 
and allows our measurements to reach their full constraining 
power. 
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